Label-free detection of small and large molecule interactions, and activities in biological systems

ABSTRACT

A system for quantitative detection and analysis of the interactions of molecules with molecular receptors on the surfaces of biological cells based on detecting a mechanical deformation in the membrane of a cell associated with the molecular interactions, which works for both large and small molecules. The mechanical deformation can be detected with high precision in real time from an optical image of the cell with a differential detection method. The system can be also used to detect the electrical activities, such as ion channel opening and closing, as well as action potential propagation in neuronal cells.

CROSS-REFERENCE TO RELATED APPLICATION

This application is a continuation of co-pending U.S. application Ser. No. 14/968,331 filed Dec. 14, 2015 having the same title which is a non-provisional application of expired U.S. Provisional Patent Application No. 62/091,828, filed Dec. 15, 2014, entitled “LABEL-FREE DETECTION OF SMALL AND LARGE MOLECULE INTERACTIONS, AND ACTIVITIES IN BIOLOGICAL SYSTEMS,” to the same inventors herein and claims the priority benefit of that filing date. U.S. application Ser. Nos. 14/968,331 and 62/091,828 are incorporated by reference.

TECHNICAL FIELD

The present invention relates to the field of biomedical research, screening and development of drugs, and discovery, validation and detection of biomarkers for diseases, and, more particularly, to a method based on nanometer-precision tracking of molecular binding-induced mechanical deformation in the cell membrane via optical imaging, which is label-free, real time and non-invasive.

BACKGROUND

Measuring molecular binding or interactions is a basic task for understanding many biological processes and for developing relevant applications. A particularly important application is to determine the binding properties of drugs with their corresponding membrane protein receptors, which are the largest type of drug targets.' Molecular binding is quantified by kinetic constants, which have, thus, become a criterion in preclinical drug screening.^(2,3)

Advances in structural biology have led to an exponential growth in the number of membrane proteins with determined 3D structures. However, as noted above, in order to understand the cellular functions of membrane proteins, it is also necessary to determine the interaction kinetics of the membrane proteins with various molecules. This is because cells perform many functions, including communication, via the interactions of their membrane proteins with molecules in the extracellular medium. A capability to quantify membrane protein interactions with molecules is also critical for discovering and validating drugs because most drug targets are membrane proteins. Despite the importance, developing such a capability that can measure the interactions of molecules with membrane proteins in the natural lipid environment has been a difficult task.

One traditional method for determining the kinetic constants is to extract molecular receptors from cells, immobilize them on a solid surface after purification, and then expose them to the drug for binding.⁴ Although useful, such methods can be problematic, especially when the receptors are membrane proteins,⁵ which currently count for more than a half of the drug targets.⁶ Due to their unique amphiphilic structures, it is difficult to ensure that the purified membrane proteins retain their native structures and functions.⁴ Because of the heterogeneous nature of cells, it is also important to study each of the individual cells. These capabilities, if developed, will benefit not only drug discovery, but also drug resistance study, which is a common but difficult problem in medicine.⁷⁻⁹

Typically, methods for studying molecular interactions use radioactive or fluorescent labels. These end point assays do not provide kinetic constants that are needed to quantify the membrane interactions and functions. To determine the kinetic information, the current practice involves extracting membrane proteins from cells, purifying them from the extracts, immobilizing the purified proteins on a solid surface, and then exposing them to a ligand for kinetic study. The procedures are not only laborious, but also prone to alteration of the native functions of membrane proteins, especially integral membrane proteins that are permanently attached to the membrane. Furthermore, the isolation of membrane proteins from their native cellular environment prevents one from studying the allosteric effect in the molecular interactions, and examining heterogeneous nature of cells. A more serious limitation of the existing technologies is that the detection signal diminishes with the mass of the molecule, making them difficult for detecting small molecules, which play many important roles in cellular functions, and represent the vast majority of the existing drugs.

Various methods have been developed for in situ measurement of drug-receptor binding. A popular method is kinetic exclusion assay, which measures the concentration of free drugs remaining in the supernatant after the binding equilibrium in cell suspension is achieved with a labeled detection technology.¹⁵⁻¹⁷ This is an end-point assay, and not suitable for extracting the kinetic constants, including the association and dissociation constants. Furthermore, the use of labels is not only labor intensive but may also affect the native binding behaviors of the molecular receptors. Label-free technologies, such as quartz crystal microbalance, have been developed for studying drug binding properties.¹⁸ Although useful, they lack spatial resolution required to study the variability between different individual cells, map heterogeneous distribution of receptors in the cell membrane, and distinguish non-specific binding onto the sensor surface from specific binding to the receptors on the cells. Another label-free detection technology is surface plasmon resonance (SPR) technique that can monitor the lectinglycoprotein interactions in single cells.¹⁹ However, like the quartz crystal microbalance, SPR signal is proportional to the mass of the molecule (e.g., drug), which has limited sensitivity for detecting drug molecules with typically small molecular masses.

The present invention overcomes the shortcomings discussed above and, for the first time, discloses a system that can detect the binding of both large and small molecules with the molecular receptors in single cells, and analyze the corresponding binding kinetic constants. The method can be applied to measure the binding of drugs with their membrane receptor targets in the native cellular membranes, and to analyze cell-to-cell variability of the binding kinetics by measuring mechanical deformation of cells upon interactions of the cellular membrane proteins with molecules in the extracellular medium. A capability of real time analysis of the interactions in single cells by analyzing the mechanical deformation with sub-nm resolution is also provided for the first time. For small molecules, the present method represents the first kinetic measurement while the equilibrium constants extracted from the present method are consistent with those obtained with endpoint radioactive labeling assay. The imaging capability allows revelation of cell-to-cell variability of difference cells, and region-to-region variability within the same cell. The detection principle of the present invention may also be used to monitor the electrical activities in neurons.

BRIEF SUMMARY OF THE DISCLOSURE

This summary is provided to introduce a selection of concepts in a simplified form that are further described below in the Detailed Description. This summary is not intended to identify key features of the claimed subject matter, nor is it intended to be used as an aid in determining the scope of the claimed subject matter.

A system for quantitative detection and analysis of the interactions of molecules with molecular receptors on the surfaces of biological cells based on detecting a mechanical deformation in the membrane of a cell associated with the molecular interactions is disclosed. An imaging apparatus captures a time sequence of images of a biological object includes an injector coupled to introduce a substance that interacts with the biological object. A processor is coupled to receive data from the imaging apparatus and adapted to determine the mechanical deformation of the biological object associated with the interaction of the substance with the biological object from the images, the processor being further adapted to analyze the interaction from the mechanical deformation.

BRIEF DESCRIPTION OF THE DRAWINGS

While the novel features of the invention are set forth with particularity in the appended claims, the invention, both as to organization and content, will be better understood and appreciated, along with other objects and features thereof, from the following detailed description taken in conjunction with the drawings, in which:

FIG. 1A schematically shows illumination of the setup for edge tracking detection.

FIG. 1B illustrates differential optical detection for accurate tracking of cell edge change induced by analyte receptor interaction.

FIG. 1C graphically illustrates a typical binding curve as determined from the cell edge movement.

FIG. 1D graphically illustrates an example where the root mean square of the fixed cell edge change is about 0.46 nm.

FIG. 1E is a cartoon that illustrates cell edge changes over time during the binding process. i, ii, and iii correspond to the stages marked in FIG. 10.

FIG. 2A shows WGA interaction with glycoproteins featuring phase contrast images of fixed CP-D cells for 20 g/ml WGA binding.

FIG. 2B shows WGA interaction with glycoproteins featuring phase contrast images of fixed CP-D cells for 5 g/ml WGA binding.

FIG. 2C graphically illustrates an example for averaged cell edge movement over the whole cell for 20 g/ml WGA.

FIG. 2D graphically illustrates an example for averaged cell edge movement over the whole cell for 5 g/ml WGA.

FIG. 3A shows acetylcholine interaction with nicotinic acetylcholine receptors in cells, specifically a phase contrast image of the fixed human a4β2 transfected SH-EP1 cells, where the white arrow marks the cell under analysis.

FIG. 3B shows averaged cell edge movement over the whole cell and fitting results during the binding process for acetylcholine of different concentrations.

FIG. 3C graphically illustrates an example wherein the equilibrium constant (KD) was determined to be 2.601 μM by plotting the edge movement vs. acetylcholine concentration and fitting the data with the Langmuir isotherm.

FIG. 3D shows acetylcholine interaction with nicotinic acetylcholine receptors in cells, specifically a phase contrast image of the fixed human a4β2 transfected SH-EP1 cells, where the white arrow marks the cell under analysis.

FIG. 3E shows averaged cell edge movement over the whole cell and fitting results during the binding process for acetylcholine of different concentrations.

FIG. 3F graphically illustrates an example wherein the equilibrium constant (K_(D)) was determined to be ˜26 nM by plotting the edge movement vs. acetylcholine concentration and fitting the data with the Langmuir isotherm.

FIG. 3G illustrates phase contrast images of wild type SH-EP1 cells.

FIG. 3H represents response of wild type SH-EP1 cells to 500 nM acetylcholine.

FIG. 4A and FIG. 4C show heterogeneity of small molecule interactions with cell membrane receptors. Phase contrast images of fixed human a4β2 transfected SH-EP1 cells, where the numbers mark the cells under analysis. Scale bar: 10.

FIG. 4B graphically illustrates an example of binding kinetics of cells 1 and 2 (50 μM Acetylcholine).

FIG. 4D graphically illustrates an example of binding kinetics of cells 3 and 4 (50 μetylcholine).

FIG. 4E graphically illustrates an example of binding kinetics at different locations (50 μetylcholine).

FIG. 4F shows a differential image of the cell shown in the insets of FIG. 4E obtained by subtracting the image recorded at 0 s from that at 68 s.

FIG. 5A and FIG. 5C show phase contrast images of fixed human a4β2 transfected SH-EP1 cells.

FIG. 5B and FIG. 5D graphically illustrate binding kinetics of cells.

FIG. 5E graphically illustrates binding kinetics at different locations of cell in f (100 nM Acetylcholine).

FIG. 5F shows a phase contrast image of fixed human a4β2 transfected SH-EP1 cell.

FIG. 6A-FIG. 6E schematically show a mechanism for trapping cells onto the holes of a microplate.

7A shows a bright field image of a floating cell clamped with glass micropipette.

FIG. 7B shows a cell modulated by an external applied AC field.

FIG. 8A schematically shows the setup of z-modulation FFT microscopy.

FIG. 8B schematically shows the application of temporal Fast Fourier Transform (FFT) applied to the bright field image sequences over a time window.

FIG. 8C schematically shows periodically modulated image intensity of a pixel, and its corresponding FFT spectrum, where a peak is located at the frequency of modulation.

FIG. 8D and FIG. 8E illustrate individual 42 nm polystyrene nanoparticles that are invisible in bright field image, but are clearly visible in the modulation image.

FIG. 8F and FIG. 8G contrast a z-modulation FFT image showing superior contrast compared with the traditional bright field image obtained with same camera.

FIG. 9A and FIG. 9B show images from a preliminary test of WGA binding kinetics on a SHEP1 cell with the z-modulation FFT microscope are shown.

FIG. 9C illustrates binding kinetic plots of 10, 20 and 50 μg/mI WGA.

In the drawings, identical reference numbers identify similar elements or components. The sizes and relative positions of elements in the drawings are not necessarily drawn to scale. For example, the shapes of various elements and angles are not drawn to scale, and some of these elements are arbitrarily enlarged and positioned to improve drawing legibility. Further, the particular shapes of the elements as drawn, are not intended to convey any information regarding the actual shape of the particular elements, and have been solely selected for ease of recognition in the drawings.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

The following disclosure describes several embodiments for metabolic analyzers that are based on detection of several metabolic signatures. Several features of methods and systems in accordance with example embodiments are set forth and described in the figures. It will be appreciated that methods and systems in accordance with other example embodiments can include additional procedures or features different than those shown in the figures. Example embodiments are described herein with respect to a portable metabolic analyzer system. However, it will be understood that these examples are for the purpose of illustrating the principles, and that the invention is not so limited. Additionally, methods and systems in accordance with several example embodiments may not include all of the features shown in the figures.

Unless the context requires otherwise, throughout the specification and claims which follow, the word “comprise” and variations thereof, such as, “comprises” and “comprising” are to be construed in an open, inclusive sense that is as “including, but not limited to.”

Reference throughout this specification to “one example” or “an example embodiment,” “one embodiment,” “an embodiment” or combinations and/or variations of these terms means that a particular feature, structure or characteristic described in connection with the embodiment is included in at least one embodiment of the present disclosure. Thus, the appearances of the phrases “in one embodiment” or “in an embodiment” in various places throughout this specification are not necessarily all referring to the same embodiment. Furthermore, the particular features, structures, or characteristics may be combined in any suitable manner in one or more embodiments.

Definitions

Generally, as used herein, the following terms have the following meanings when used within the context of sample collection or analysis:

As used herein, “plurality” is understood to mean more than one. For example, a plurality refers to at least 3, 4, 5, 70, or more.

As used herein, “neuron” refers to a nerve cell that receives and sends electrical signals over long distances within the body. A neuron receives electrical input signals from sensory cells (called sensory neurons) and from other neurons. The neuron sends electrical output signals to muscle neurons (called motoneurons or motor neurons) and to other neurons.

As used herein, “MADMI” refers to a method for mechanically amplified detection of molecular interactions that provides amplified molecular interaction signals, allowing measurement of the binding kinetics of molecules with membrane proteins on single live cells, and analysis of the heterogeneous nature of the binding kinetics between different cells, and different regions of a single cell.

As used herein, “Nicotinic acetylcholine receptors,” or “nAChRs,” are neuron receptor proteins that signal for muscular contraction upon a chemical stimulus.

Principle of the invention

The basic detection principle of the present invention is to measure a deformation (shape or size) in the cell that accompanies the binding of molecules with the molecular receptors on the cell surface. This deformation occurs because the law of thermodynamics predicts that when molecules bind to the surface the surface tension changes. According to thermodynamics, the relative surface concentration of molecules bound on the membrane surface is given by

$\begin{matrix} {{\Gamma = {- \left( \frac{\delta\gamma}{\delta\mu} \right)_{T,P,}}},} & (1) \end{matrix}$

where γ is the surface tension and μ is the chemical potential of the molecules. For ideal solutions, the chemical potential is related to the bulk concentration, c, according to

dμ=k _(B) TdΓ(Inc),   (2)

where k_(B) is the Boltzmann constant, and Tis temperature. Substituting Eq. 2 into Eq. 1, we have

dγ=kBTdΓ(dc/c)   (3)

If we assume the molecular binding follows the Langmuir isotherm, Eq. 3 can be simplified as

$\begin{matrix} {{{d\; \gamma} = {{RT}\left( \frac{\Gamma_{0}}{\Gamma_{0} - \Gamma} \right)}},} & (4) \end{matrix}$

where Γ₀ is the surface concentration at the full (maximum) coverage. When Γ<<Γ₀, Eq. 4 can be further simplified and takes the form of,

dγ=k_(B)TdΓ.   (5)

Eq. 5 clearly shows that the surface tension of the cell membrane changes with the molecular binding.

Some of the receptor proteins are ion channels, which open and close to allow different ions to move in and out of a cell. Associated with the ion re-distribution is an electrical potential change across the cell's membrane, particularly neurons, giving rise to electrical activity, such as action potential that is critical for information processing in brains. The change in the ion distribution or electrical activities also leads to a mechanical deformation in the cell, as described the Lippmann equation, given by

$\begin{matrix} {{q = \left( \frac{\delta\gamma}{\delta \; V} \right)_{T,p,\mu}},} & (6) \end{matrix}$

where q is surface charge density and V is the electrical potential. Eq. 6 predicts that for a given charge density, a change in the electrical potential results in a change in the surface tension, thus a mechanical deformation in the cell. This is the principle for the detection of electrical activities of neurons, such as action potential, by measuring and analyzing the cell deformation.

Note that Eq. 5 was derived with many assumptions, but the basic concept that a molecular binding event causes a change in the surface tension is expected to be a general phenomenon. However, the actual amount of conformational change in the cell depends on the biding strengths, as well as the geometry and mechanical properties of the cells. For practical applications, one typically wants to determine the binding kinetic constants, which can be obtained from the time dependence of the cell membrane deformation (mechanical response) during the association (binding) and dissociation (unbinding) processes. In other words, there is no need to relate the measured mechanical response to the microscopic scaled processes.

Description of the System

Referring now jointly to FIG. 1A-FIG.1E, an overview of edge tracking detection employing the methods disclosed herein is illustrated. To detect a small amount of molecules or a weak molecular binding event, it is critical to be able to measure a small amount of mechanical deformation in the cell. This task is accomplished in the present system by detecting and tracking the edge movement of a cell image.

Referring now particularly to FIG. 1A, there shown schematically is an example of an illumination setup. The setup includes a cell 5 with receptors 15 on a glass slide 22. Analytes 20 are introduced and a differential interference contrast (DIC) microscope 10 is positioned to track activity. A traditional bright field optical microscopy may be used, but better image contrast can be obtained with phase contrast microscopy, differential interference contrast (DIC) microscopy, and various interference optical microscopic techniques. In one useful example an inverted phase contrast microscope with 40× phase 2 objective 10 was used for detection. The DIC microscope 10 transmits time sequenced images to a sensor 41, which in turn passes data representing the time sequenced images to a processor 43 for analysis. The sensor 41 may be any commercially available device such as a CCD video camera or equivalent. The processor may advantageously include a CPU or equivalent device for processing digital data, for example. The data for time sequenced images may be analyzed to determine differences in the edge boundaries for the cell at different times in response to analytes. Edge boundary measurement with software programs is well known. Using the microscope in combination with the sensor and processor detecting a mechanical deformation in the membrane of a cell associated with the molecular interactions can be performed as desired.

Referring now particularly to FIG. 1B, differential optical detection for accurate tracking of cell edge change induced by analyte receptor interaction is illustrated. The edge of a cell 5 can be manually defined from the optical image of the cell, but an imaging processing algorithm can automatically detect the cell edge, providing high throughput. After determining the cell edge, a region of interest (ROI) is selected at a location of the cell edge. The ROI can be a rectangular shape30 with the cell edge 32 passing through the center portion of the rectangle, dividing the rectangle into two equal halves, one half 34 is inside of the cell, and the second half 36 falls outside of the cell. We denote the intensities of the two halves of the differential image as A and B, and (A−B)/(A+B) is calculated and used to determine the movement of the cell edge at each location. We refer (A−B)/(A+B) to as differential image intensity, and this method of edge movement tracking differential optical detection.

An alternative approach is to first determine a differential image within the ROI by shifting the image within the ROI perpendicular to the cell edge by a small distance (e.g., 1 pixel), and then calculating the difference between the shifted image and the original image pixel by pixel. We denote the intensities of the two halves of the differential image as dA and dB, and (dA−dB)/(dA+dB) is calculated and used to determine the movement of the cell edge at a given location. This method of edge movement tracking is also referred to as differential optical detection.

The relation between the cell edge movement and the measured differential image intensity change from one of the above differential optical detection method can be determined as below. The pixel density of each image is enlarged 5 times by adding additional pixels with a bilinear interpolation approach. The distance between two pixels in the interpolated image corresponds to a real distance of the object (cell) is 37 nm (Pike F032B CCD, Allied Vision Technologies, Stadtroda, Germany). The edge of one cell is manually chosen and the centroid (O) of the cell is determined (FIGS. 1B and 1C). A polar coordinate system is setup with the centroid serving as the pole. The cell edge movement was calculated at every 10 starting from 0 (FIG. 1E). The ROI at a certain point of the cell edge (point A) is then shifted by different numbers of pixels outwards (perpendicular to the tangential line at point A), and the corresponding changes in the differential image intensity is determined from the image. The relation between the differential image intensity and the cell edge movement (pixels) is found to be linear within a certain range, which serves as a calibration curve to determine the cell edge movement (mechanical deformation) from the differential image intensity. This differential optical detection method can accurately detect cell edge movement with a detection limit as small as 0.5 nm, the size of an atom (FIG. 1D).

EXAMPLES

Referring now jointly to FIG. 2A-FIG. 2D, there shown are illustrative drawings of WGA interaction with glycoproteins. White arrows 40, 42 mark the cells under analysis in FIG. 2A and FIG. 2B respectively. The scale bar for FIGS. 2C and 2D is 10 μm. FIG. 2C graphically illustrates an example for averaged cell edge movement over the whole cell for 20 g/ml WGA. FIG. 2D graphically illustrates an example for averaged cell edge movement over the whole cell for 5 g/ml WGA. The scale bar for FIGS. 2C and 2D is 10 μm. The Y axis for FIGS. 2C and 2D represents edge movement in nm ranging from −5 to 25 nm. The X axis represents time in seconds.

The following examples demonstrate the detection of the binding of both large and small molecules with membrane receptors in cells using the system and method in this system. An inverted microscope (Olympus X81) equipped with phase 2 condenser and phase 2 40× objective was used with illumination from the top of the sample cells.

Large Molecule-Membrane Receptor Protein Interactions

To demonstrate the capability of the system and method for detecting and analyzing molecular binding of membrane receptors of cells, the binding kinetics of wheat germ agglutinin (WGA) and glycoprotein on Barrett's esophagus derived CP-D (CP-18821) cells was studied. WGA is one kind of lectins that binds to N-Acetyl glucosamine (GlcNAc) and sialic acid groups. The CPD cells were fully attached onto a glass slide.

FIGS. 2A and 2B show the phase contrast images of attached fixed CP-D cells. The measurement was carried out by flowing 1× PBS over the cells for 30 s with a flow rate of 350 L/min to obtain the baseline. At time 0 s, WGA solution in 1× PBS was introduced for 90 s allowing the association (binding) of WGA with CP-D glycoprotein on the cell surface. During the association process, the cell edge moves outwards as shown in FIGS. 2C and 2D. Then WGA solution was switched to 1× PBS to allow the bound WGA to dissociate from the CP-D cells at time 90 s. FIGS. 2C and 2D show that the cell edge moves back to the original position during the dissociation process. By fitting the data with the first order kinetics, association rate constants (k_(on)), dissociation rate constants (k_(off)), and dissociation constant (K_(D)) were found to be k_(on)=1.4×10⁵ M⁻¹ s⁻¹, k_(off)=2.9×10⁻³s⁻¹ s⁻¹, K_(D)=0.021 μM for 20 ug/ml WGA and k_(on)=1.09×10⁵ M⁻¹ s⁻¹, k_(off)=1.5 6×10⁻³ s⁻¹, K_(D)=0.014 μM for 5 ug/ml WGA.

Small Molecule-Membrane Receptor Protein Interactions

In order to demonstrate small molecule binding to single cells with the present system, the binding of acetylcholine with nicotinic acetylcholine receptors (nAChRs) was studied. Engineered SH-EP1 cells that expressed human a4β2 receptors were used to examine the binding kinetics of acetycholine with nAChRs. Human a4β2 mainly exists in the brain and is related to the nicotine addiction. Study of the nAChRs and small molecule interactions is important for revealing the mechanism of nicotine addiction as well as development of drugs to treat nicotine addiction.

Referring now jointly to FIG. 3A-FIG. 3C, acetylcholine interaction with nicotinic acetylcholine receptors in cells, specifically a phase contrast image of the fixed human a4β2 transfected SH-EP1 cells. The phase contrast image of fixed SH-EP1-ha482 cells is shown in FIG. 3A, where the white arrow 50 marks the cell under analysis. The scale bar is 10 μm.

1× PBS buffer was continuously flowing over the cell for 25 seconds to obtain a stable baseline. At time 0 s, the buffer was changed to acetylcholine solution of different concentrations in 1× PBS. After association, the acetylcholine solution was switched back to 1 × PBS to allow dissociation.

Referring now particularly to FIG. 3B averaged cell edge movement over the whole cell and fitting results during the binding process for acetylcholine of different concentrations is shown. The Y axis represents edge movement in nm ranging from about 0 to 40 nm. The X axis represents time in seconds. Averaged cell edge movement over the whole cell is indicated by the experimental data points 70 fitted results are shown as a plurality of curves 72. From bottom to top the buffer (negative control), was 0.5 μM, 2 μM, 5 μM, and 50 μM.

As shown in FIG. 3B, the cell edge moves outwards during the association phase and retracts back during the dissociation phase. It also shows that the amount of cell edge movement during the association process increases with the acetylcholine concentration, which is expected for first-order reaction kinetics. The association rate constant (k_(on)) and dissociation rate constant (k_(off)) were found to be 9.4×10 ³ M⁻¹ s^(−1A)nd 2.47×10⁻² s⁻¹ respectively, from which the dissociation constant (K_(D)) was determined to be 2.7 μM.

Referring now particularly to FIG. 3C, there illustrated is an example wherein the equilibrium constant (K_(D)) was determined to be 2.601 μM by plotting the edge movement vs. acetylcholine concentration and fitting the data with the Langmuir isotherm. The Y axis represents edge movement in nm ranging from about 0 to 40 nm. The X axis represents concentration (M) ranging from about 10⁻⁷ to 10⁻⁴ M.

Referring now jointly to FIG. 3D-FIG. 3F, acetylcholine interaction with nicotinic acetylcholine receptors in cells, specifically a phase contrast image of the fixed human a4β2 transfected SH-EP1 cells. The phase contrast image of fixed SH-EP1-α4β2 cells is shown in FIG. 3D, where the white arrow 350 marks the cell under analysis. The scale bar is 20 μm. 1× PBS buffer was first introduced to flow over the cell for 25 s, and then the buffer was switched to an acetylcholine solution in 1× PBS. After association, the acetylcholine solution at each concentration was switched back to 1× PBS to allow for dissociation. The above procedure was repeated for different acetylcholine concentrations.

As shown in FIG. 3E, the cell edge expands during the association phase and retracts during the dissociation phase. FIG. 3E also shows that the amount of cell expansion during the association process increases with the acetylcholine concentration, which is expected for first-order binding kinetics. The association (kon) and dissociation (koff) rate constants were found to be 1.2×10⁻⁶ M⁻¹s⁻¹ and 2.2×10⁻² s⁻¹, respectively, which represent the first direct measurement of the kinetic constants for the binding of the neurotransmitter to the nAChRs in intact cells.

From kon and koff, the equilibrium dissociation constant (KD=koff/kon) was determined to be 18.1 nM. By plotting the equilibrium response versus acetylcholine concentrations (FIG. 3F), the equilibrium constant (KD) was found to be ˜26 nM, which is consistent with that obtained by kinetics measurement. As this is the first kinetic measurement of acetycholine binding to nAChRs, the findings cannot be compared to other reference technologies or prior data. However, the equilibrium dissociation constant determined here is in good agreement with the average Ki determined with radioligand binding assay, which involved centrifuge and formation of cell pellets.

Referring now to FIG. 3G, phase contrast images of wild type SH-EP1 cells G,

H are shown. The scale bar is 20 μm. As a control experiment, the measurement was carried out with wild type SH-EP1 cells, which do not have nAChRs expressed on the cell surfaces, and observed no deformation in the cell membrane. This result demonstrated that the mechanical deformation in the engineered SH-EP1 cells was indeed due to the specific binding of acetylcholine to the expressed nAChRs.

Referring now to FIG. 3H, response of wild type SH-EP1 cells to 500 nM acetylcholine is graphically shown. Arrow 350 marks the switch from 1× PBS to 500 nM acetylcholine and the arrow 352 indicates the change of solution back to 1× PBS. Curve 370 represents the measurements for cell G and lighter curve 370 represents the measurements for cell H.

Heterogeneity

Referring now jointly to FIG. 4A-FIG. 4F, heterogeneity of small molecule interactions with cell membranes are illustrated. The present system makes it possible to examine the cell-cell variation and heterogeneity within a cell. FIGS. 4A and 4C are the phase contrast images of fixed SH-EP1-hα4β2 cells and numbers in circles mark the cells under analysis. The scale bar is 10 μm. Cell 1 and cell 2 are on the same glass slide, and cell 3 and cell 4 are on a different glass slide. The responses of four cells to 50 -μM acetylcholine are shown in FIGS. 4B and 4D, and the corresponding kinetic constants are given in Table 1, which show significant differences in the binding kinetics. Different regions on a cell edge also show variations in the binding kinetics (FIGS. 4E and 4F). The scale bar is 10 μm. Insets 55, 65 show phase contrast images of fixed human α4β2 transfected SH-EP1 cells recorded at the time 0 s and 68 s, respectively.

Referring particularly to FIG. 4B, binding kinetics of cells 1 and 2 are graphed. The Y axis represents cell edge movement in nm and the X-axis represents time in seconds. Curve 470 represents fitted results for experimental data 472 for cell 1. Curve 473 represents fitted results for experimental data 474 for cell 2. Similarly, FIG. 4D shows binding kinetics of cells 3 and 4. The Y axis represents cell edge movement in nm and the X-axis represents time in seconds. Curve 483 represents fitted results for experimental data 484 for cell 3. Curve 480 represents fitted results for experimental data 482 for cell 4.

TABLE 1 Association rate constants (k_(on)), dissociation rate constants ₍k_(off),) and equilibrium constants (K_(D)) for four cells as shown in FIGS. 4A-4D. k_(on) (M⁻¹s⁻¹) k_(off) (s⁻¹) K_(D)(μM) Cell 1 1.00 × 10³ 0.0152 15.1 Cell 2 5.02 × 10² 0.0139 28.1 Cell 3 4.01 × 10² 0.0142 36.1 Cell 4 1.26 × 10³ 0.0150 11.9

Referring now to FIG. 5A and FIG. 5C phase contrast images of fixed human α4β2 transfected SH-EP1 cells are shown. Cells under analysis include cells 5-7. The scale bars for both figures is 20 μm.

Referring to FIG. 5B, binding kinetics of cells 5 and 6 are graphed. The Y axis represents cell edge movement in nm and the X-axis represents time in seconds. Curve 570 represents fitted results for experimental data 572 for cell 5. Curve 573 represents fitted results for experimental data 574 for cell 6. Similarly, FIG. 5D shows binding kinetics of cells 7 and 8. The Y axis represents cell edge movement in nm and the X-axis represents time in seconds. Curve 583 represents fitted results for experimental data 584 for cell 7. Curve 580 represents fitted results for experimental data 582 for cell 8.

TABLE 2 Association rate constants (kon), dissociation rate constants (koff), and equilibrium constants (KD) for four cells as shown in FIG. 5A-FIG. 5D. kon (M⁻¹s⁻¹) koff (s⁻¹) KD (nM) Cell 5 7.32 × 10⁵ 0.0161 22.0 Cell 6 8.33 × 10⁵ 0.0217 26.1 Cell 7 3.42 × 10⁵ 0.0298 87.1 Cell 8 5.76 × 10⁵ 0.0172 29.9

Referring to FIG. 5E binding kinetics at different locations of cell (100 nM Acetylcholine) is graphically illustrated. The Y axis represents cell edge movement in nm and the X-axis represents time in seconds.

Referring to FIG. 5F, a phase contrast image of fixed human α4β2 transfected SH-EP1 cell is shown. The amount of cell membrane deformation after acetylcholine binding ranges from about 80 nm, for example, at arrow 590 to about 9−40 nm, at, for example, region 592 (scale bar: 20 μm).

Experimental Details Materials

Wheat germ agglutinin (WGA) and acetylcholine chloride were purchased from Sigma-Aldrich (St. Louis, Mo.). 1× phosphate buffer saline (PBS) pH=7.4 was used as buffer for all binding experiments. All samples were prepared in 1× PBS buffer.

Cell Culture

The human α4β2 transfected human epithelial SH-EP1 cells were cultured in a humidity incubator at 37 with 5% CO2 and 70% relative humidity. Dubelco's Modified Eagle's Medium (DMEM, Lonza, Walkersville, Md.) with 10% Fetal Bovine Serum (FBS, Life Technologies, Carlsbad, Calif.) and penicillin-streptomycin (BioWhittaker, Basel, Switzerland) were used as culture medium. SH-EP1-h α4β2 cells were cultured in 25 cm² flask until approximately 80% confluence was reached for passage. 0.05% trypsin-EDTA (Life Technologies, Carlsbad, Calif.) was used for cell passage.

The CP-D cell were cultured in an incubator at 37 with 5% CO2 and 70% relative humidity. Cells were cultured in 25 cm² flask with 1× Keratinocyte-SFM (Life Technologies, Carlsbad, Calif.) and penicillin-streptomycin (BioWhittaker, Basel, Switzerland) as culture medium. When cells reached approximately 80% confluent, cells were passaged with 0.05% trypsin-EDTA (Life Technologies, Carlsbad, Calif.).

For experiments, cells were cultured overnight on the bare glass slides (22×60 mm micro cover glass, VWR, Radnor, Pa.) in a silicone well (FlexiPERM, Greiner bio-one, Monroe, N.C.) placed on top of it in order to let cells attach on the surface. Cells on glass slides were also cultured in the incubator at 37 with 5% CO2 and 70% relative humidity. Cells were incubated in 4% paraformaldehyde for 10 μmin at room temperature for fixation and then ready for experiments. Before the measurement, the small silicone well was changed to a homemade PDMS well with 2 cm in length, 1 cm in width and 1 cm in height.

High Throughput Analysis of Individual Cells with Minimized Micro-Motions

As discussed in detail above, thermodynamic principles (Eqs. 1 and 2) predict that a mechanical deformation always accompanies a molecular binding on a cell surface. This is the basic principle of MADMI. However, thermodynamics does not bring to light how much a cell will deform for a given molecular binding. This information is not needed for binding kinetic constant measurement because the kinetics constants can be determined from fitting the relative mechanical deformation vs. time plot with a kinetics model. Precise prediction of the amount of mechanical deformation associated with a binding event would require molecular scaled knowledge about the nature of the molecule-membrane protein interactions, as well as the surface density of the membrane proteins, and mechanical property of the cell. Obtaining such detailed knowledge and analyzing MADMI data in terms of the microscopic knowledge are beyond the scope of the present project. However, it is believed that the amount of mechanical deformation in MADMI increase with the binding strength and surface density of the target membrane proteins.

Referring now to FIG. 6A, trapping of cells onto the holes of a microplate is schematically shown. To minimize micro-motion of live cells and achieve high throughput analysis of molecular binding kinetics, a method to trap individual cells 601 for real-time study of molecular binding kinetics of membrane proteins using MADMI will use a microplate 600 with an array of μm-scaled holes 605, such that individual cells 601 can be trapped onto the holes with an applied pressure difference across the holes. Preliminary results show that MADMI works well with adherent cells. At least one cell 601 will be trapped onto a hole 603 via a negative pressure created with a laminar flow 605 underneath the microplate 600. An objective lens 612 may be used for obtaining images. A similar method has been used in commercial automated patch clamping systems, where micron-sized holes are created in the bottom of the wells of a microplate.

Referring now jointly to FIG. 6B-FIG. 6E, apparatus for trapping cells is illustrated. FIG. 6B is a photo of a 96 well automated patch clamp plate 600 of an lonFlux Automated Patch Clamp System. FIG. 6C is an expanded view showing the fluidic design and function of each well in the lonFlux Automated Patch Clamp System of FIG. 6B. FIG. 6D is a more expanded view of a portion of the lonFlux Automated Patch Clamp System showing the cell-trapping zone 610. FIG. 6E shows an image of trapped cells. lonFlux Automated Patch Clamp Systems are available from Fluxion Biosciences, Inc., South San Francisco, Calif. 94080.

Referring now to FIG. 7A, a bright field image of a floating cell clamped with glass micropipette is shown. A cell 701 is physically connected with micropipette 704 in a whole-cell patch configuration. The white spot 705 in the center of the cell is the tip of the micropipette. As a preliminary test to examine the feasibility of the method, a cell was trapped with a micropipette used in patch clamp. An AC electrical field having a frequency of 37 Hz was then applied across the cellular membrane to create a periodic mechanical deformation in the membrane.

Referring now to FIG. 7B, a cell 701 modulated by an external applied AC field from micropipette 704 is shown. The electromechanical coupling induced mechanical deformation was precisely measured by MADMI method and visualized with nanometer precision ranging from 0 to about 5 nm. The mechanical deformation occurred because of an electromechanical coupling effect, which was precisely measured by MADMI with nanometer precision. This preliminary test shows that it is possible to trap a cell to reduce its micro-motions, and to detect the mechanical deformation of the cell membrane with nanometer precision.

With the trapping of the cells, cellular micro-motions may be deteremined by performing noise spectrum analysis, and studying the micro-motions by applying different negative pressure, and effect of the size of the holes. MADMI may then be applied to determine molecular interaction kinetics using both large and small molecules, and compare the results with those of the adherent cells. Based on the finding, the cell trapping method may be optimized by tuning the diameter, flow rate and shear force for rapid cell trapping, in order to minimize the micro-motions and maximizing the throughput.

A low-noise imaging technology for studying low-density membrane proteins with MADMI. The density or abundance of the membrane proteins on cell surface varies over several orders of magnitudes (˜1 0³-1 0⁸ protein/cell, or 1-10⁵ protein/μm²). Since the binding signal is proportional to the membrane protein receptor density, in order to study low-density membrane proteins with MADMI, it is necessary to improve the detection limit of the cell membrane deformation. To achieve this goal, a low noise z-modulation Fast Fourier Transform (FFT) optical imaging technology may be used as described below. The FFT technology will periodically modulate the sample along the optical axis (z-axis) to create a periodic modulation in the focus, and determines the associated image contrast modulation with FFT, which minimizes noise and provides superior detection limit. We expect to reach a detection limit of 0.25 nm for the membrane deformation. Such a detection limit will help detect the bindings of molecules to low-density membrane proteins, e.g., 10 proteins/μm².

The z-modulation FFT microscope will improve the signal to noise ratio by modulating the sample stage vertically (z-axis). This is possible because the z-modulation FFT microscope will detect signals associated with the sample only, and it reduces noise from light source, optics and camera, and also removes effects due to defects or dirt in the optics of the microscope. Additionally, the periodic modulation will allow effective use of FFT filter to remove all types of noise that has different frequencies from the modulation frequency. The z-modulation microscope may advantageously be built on a traditional bright field optical microscope with an attached z-piezo modulation sample holder and proper imaging-processing algorithm. A similar method was previous demonstrated by Gineste et al. ²⁵ An important difference is that the system disclosed herein will use periodic modulation and Fourier transform algorithm to reduce noise in the image.

The basic principle of z-modulation FFT microscopy is described by the Transport of the Intensity Equation (TIE) according to the formula

$\begin{matrix} {{{{- k}\frac{\partial{I\left( {\overset{\rightarrow}{r_{\bot}},0} \right)}}{\partial z}} = {\nabla_{\bot}{\cdot \left\lbrack {{I\left( {\overset{\rightarrow}{r_{\bot}},0} \right)}{\nabla_{\bot}{\varphi \left( {\overset{\rightarrow}{r_{\bot}},0} \right)}}} \right\rbrack}}},{k = \frac{2\pi}{\lambda}}} & (3) \end{matrix}$

where I({right arrow over (r_(⊥))}, z) and ϕ({right arrow over (r_(⊥))}, z) are the intensity and phase at the position at ({right arrow over (r_(⊥))}, z) respectively, z denotes position along the optical axis and {right arrow over (r_(⊥))} denotes position within a plane normal to the optical axis, and λ is the wavelength and k is the wavenumber. This equation relates the rate of change of intensity in the direction of the optical axis to the intensity and phase of light in a plane perpendicular to the optical axis. When I({right arrow over (r_(⊥))}, 0) is constant, I₀, Eq. 3 can be simplified as

$\begin{matrix} {{\frac{\partial{I\left( {\overset{\rightarrow}{r_{\bot}},0} \right)}}{\partial z} = {{- \frac{I_{0}}{k}}\left( {\nabla_{\bot}^{2}{\varphi \left( {\overset{\rightarrow}{r_{\bot}},0} \right)}} \right)}},} & (4) \end{matrix}$

where the left of the equation is the z-modulation FFT microscopy image contrast, and the right side is the spatial Laplacian, ∇_(⊥) ², (second derivative) of image phase, ϕ({right arrow over (r_(⊥))}, 0), at the focus. Eq. 4 shows that the z-modulation FFT microscopy image contrast is proportional to the second derivative of the phase of the object. For this reason, z-modulation FFT microscopy provides excellent image contrast, which allows accurate detection of the cell deformation.

Referring now to FIG. 8A the setup of z-modulation FFT microscopy is shown. The sample, and thus the focal plane, is modulated vertically with a piezoelectric transducer driven by a periodic voltage 815, and bright field images are recorded sequentially.

Referring now to FIG. 8B, to obtain the focus modulation image, temporal Fast Fourier Transform (FFT) was applied to the bright field image sequences 820 over a time window 822. Amplitude data 830 and phase data 832 can thus be obtained.

Referring now to FIG. 8C, a periodically modulated image intensity of a pixel, and its corresponding FFT spectrum, where a peak is located at the frequency of modulation are shown. Plot 802 is a time plot of modulated image intensity. FFT intensity spectrum for a background region 804 and a sample region 805 with focus modulation at 37 Hz is shown. A peak 810 appears in the FFT spectrum at the modulation frequency for the sample region, but not in the background region. The amplitude of the peak 810 was extracted for each pixel to create a z-modulation FFT microscopy image.

Referring now jointly to FIG. 8D and FIG. 8E, to evaluate noise reduction capability of z-modulation FFT microscopy, 42 nm polystyrene nanoparticle (PSNP) were imaged on a glass surface. Individual 42 nm polystyrene nanoparticles are invisible in bright field image (FIG. 8D), but are clearly visible in the modulation image as marked by arrows 830 (FIG. 8E). Similarly, background noise in the bright field image (FIG. 8D) of cells is removed (FIG. 8G) and detailed cellular structures are enhanced in the z-modulation FFT image. FIG. 8D shows the traditional bright field image of the same sample, which cannot resolve the nanoparticles due to noise. In contrast, z-modulation FFT microscopy can clearly resolve individual 42 nm PSNPs, as shown in FIG. 8E. Another example is cultured cells. The z-modulation FFT image (FIG. 8G) shows superior contrast compared with the traditional bright field image obtained with same camera (FIG. 8F). Individual 42 nm polystyrene nanoparticles are invisible in bright field image (FIG. 8D), but are clearly visible in the modulation image as marked by arrows 830 (FIG. 8E). Similarly, background noise in the bright field image (FIG. 8D) of cells is removed (FIG. 8G) and detailed cellular structures are enhanced in the z-modulation FFT image.

Referring now jointly to FIG. 9A and FIG. 9B, images from a preliminary test of WGA binding kinetics on a SHEP1 cell with the z-modulation FFT microscope are shown. FIG. 9. WGA binding kinetic measured by z-modulation FFT microscopy. FIG. 9A Bright field image of fixed SHEP1 cell; FIG. 9B z-modulation FFT image of the same cell; FIG. 9C Binding kinetic plots of 10, 20 and 50 μg/mI WGA (from bottom to top curves). The plots were obtained by tracking the intensity changes along the edge of the modulation amplitude image. Measured kinetic constants: k_(on)=1.4×10⁴ M⁻¹s⁻¹, k_(off)=0.013 s⁻¹, K_(D)=125 nM. Scale bar: 10 μm. The kinetic curves were obtained from the amplitude of the image sequences along the edge of a cell. Kinetic constants and binding affinity were obtained by fitting the kinetic curves with first order kinetic model.

The z-modulation FFT imaging technology together with the differential optical algorithm may be used to study molecular binding kinetics of membrane proteins with different expression levels. One of the model systems planned for study is human epidermal growth factor receptor (EGFR, or HER1), which interacts with target protein drug, Panitumumab (Vectibix, FDA approved recombinant human IgG2k mAb drug manufactured by Amgen). EGFR plays an essential role in regulating normal cell signaling, and the mutation of EGFR leads to cell proliferation, angiogenesis, invasion, metastasis and inhibition of apoptosis, accounting for the pathogenesis and progression of cancer cells (33-37). Monoclonal antibodies targeting the extracellular do-main of EGFR have been used in various stages of pre-clinical development, and have shown good therapeutic efficacy for treatment of a number of cancers that have up-regulated EGFR expression level. The kinetic constants of the binding of these antibody drugs to EGFR receptors are the key parameters to characterize the efficacy of these drugs. We will use cell lines with different EGFR expression levels to measure the antibody drug interaction kinetics, and to evaluate the detection limit of MADMI.

Cell lines with different expression levels of EGFR, including A431 (high), Hela (medium) and A549 (low), and HEK293 (negative control) will be used. The EGFR expression levels of these cells have been confirmed in our lab using quantitative SPR imaging and immunofluorescence imaging. The membrane densities of EGFR in A431, Hela and A549 cells calculated from maximum SPR binding signals are 640, 270 and 140 receptors/ μm², respectively.

Binding kinetic curves of different doses of Panitumumab antibody with these four cell lines will be measured by MADMI. Kinetic constants will be obtained through globe fitting of the binding kinetics curves, and EGFR expression level will be quantified from the maximum binding signal.

The invention has been described herein in considerable detail in order to comply with the Patent Statutes and to provide those skilled in the art with the information needed to apply the novel principles of the present invention, and to construct and use such exemplary and specialized components as are required. However, it is to be understood that the invention may be carried out by different equipment, and devices, and that various modifications, both as to the equipment details and operating procedures, may be accomplished without departing from the true spirit and scope of the present invention.

REFERENCES

-   The following publications are incorporated by reference. -   1. Ehrlich, P. Address in pathology on chemotherapeutics :     Scientific principles, methods, and results. 2, 445-451 (1913). -   2. Copeland, R. A., Pompliano, D. L. & Meek, T. D.     Opinion—Drug-target residence time and its implications for lead     optimization. Nature Rev. Drug Disc. 5, 730-739 (2006). -   3. Swinney, D. C. The role of binding kinetics in therapeutically     useful drug action. Curr. Opin. Drug Disc. Develop. 12, 31-39     (2009). -   4. Ernst, R. E., High, K. N., Glass, T. R. & Zhao, Q. in Therapeutic     Monoclonal Antibodies: From Bench to Clinic (ed Zhiqiang An) (John     Wiley & Sons, Inc., 2009). -   5. Butcher, E. C. Can cell systems biology rescue drug discovery?     Nature Rev. Drug Disc. 4, 461-467 (2005). -   6. Hopkins, A. L. & Groom, C. R. The druggable genome. Nature Rev.     Drug Disc. 1, 727-730 (2002). -   7. Nahta, R. & Esteva, F. J. Herceptin: mechanisms of action and     resistance. Cancer Lett. 232, 123-138 (2006). -   8. Nahta, R., Yu, D. H., Hung, M. C., Hortobagyi, G. N. &     Esteva, F. J. Mechanisms of disease: understanding resistance to     HER2-targeted therapy in human breast cancer. Nat. Clin. Pract.     Oncol. 3, 269-280 (2006). -   9. Valabrega, G., Montemurro, F. & Aglietta, M. Trastuzumab:     mechanism of action, resistance and future perspectives in     HER2-overexpressing breast cancer. Annal. Oncol. 18, 977-984 (2007). -   10. Nagy, P. et al. Decreased accessibility and lack of activation     of ErbB2 in JIMT-1, a herceptin- resistant, MUC4-Expressing breast     cancer cell line. Cancer Res. 65, 473-482 (2005). -   11. Price-Schiavi, S. A. et al. Rat MUC4 (sialomucin complex)     reduces binding of anti-ErbB2 antibodies to tumor cell surfaces, a     potential mechanism for Herceptin resistance. Int. J. Cancer 99,     783-791 (2002). -   12. Kute, T. et al. Development of Herceptin resistance in breast     cancer cells. Cyto. Part A 57A, 86-93 (2004). -   13. Lu, Y. H., Zi, X. L., Zhao, Y. H., Mascarenhas, D. & Pollak, M.     Insulin-like growth factor-I receptor signaling and resistance to     trastuzumab (Herceptin). J. Natl. Cancer Inst. 93, 1852-1857 (2001). -   14. Nahta, R., Takahashi, T., Ueno, N. T., Hung, M. C. &     Esteva, F. J. P 27(kip1) down-regulation is associated with     trastuzumab resistance in breast cancer cells. Cancer Res. 64,     3981-3986 (2004). -   15. Darling, R. J. & Brault, P. A. Kinetic exclusion assay     technology: Characterization of molecular interactions. Assay Drug     Develop. Technol. 2, 647-657 (2004). -   16. Xie, L. et al. Measurement of the functional affinity constant     of a monoclonal antibody for cell surface receptors using kinetic     exclusion fluorescence immunoassay. J. Immuno. Methods 304, 1-14     (2005). -   17. Rathanaswami, P., Babcook, J. & Gallo, M. High-affinity binding     measurements of antibodies to cell-surface-expressed antigens. Anal.     Biochem. 373, 52-60 (2008). -   18. Pei, Z. C., Saint-Guirons, J., Kack, C., lngemarsson, B. &     Aastrup, T. Real-time analysis of the carbohydrates on cell surfaces     using a QCM biosensor: a lectin-based approach. Biosen. Bioelectron.     35, 200-205 (2012). -   19. Wang, W. et al. Label-free measuring and mapping of binding     kinetics of membrane proteins in single living cells. Nature Chem.     4, 846-853 (2012). -   20. Troise, F., Cafaro, V., Giancola, C., D'Alessio, G. & De     Lorenzo, C. Differential binding of human immunoagents and Herceptin     to the ErbB2 receptor. FEBS J. 275, 4967-4979 (2008). -   21. Pan, A. C., Borhani, D. W., Dror, R. O. & Shaw, D. E. Molecular     determinants of drug—receptor binding kinetics. Drug Disc. Today in     press (2013). -   22. DeFazio-Eli, L. et al. Quantitative assays for the measurement     of HER1-HER2 heterodimerization and phosphorylation in cell lines     and breast tumors: applications for diagnostics and targeted drug     mechanism of action. Breast Cancer Res. 13 (2011). -   23. Pearlberg, J. et al. Screens using RNAi and cDNA expression as     surrogates for genetics in mammalian tissue culture cells. Cold     Spring Harb. Symp. Quant. Biol. 70, 449-459 (2005). -   24. Hathaway, H. J. et al. Detection of breast cancer cells using     targeted magnetic nanoparticles and ultra-sensitive magnetic field     sensors. Breast Cancer Res. 13 (2011). -   25. Gineste J M, Macko P, Patterson E A, Whelan M P.     Three-dimensional automated nanoparticle tracking using Mie     scattering in an optical microscope. J Microsc-Oxford.     2011;243(2):172-8. doi: DOI 10.1111/j.1365-2818.2011.03491.x. PubMed     PMID: WOS:000292845200006. 

What is claimed is:
 1. A system for quantitative detection and analysis of the binding of molecules with molecular receptors on the surfaces of a biological object, the system comprising: a microscope positioned to image the biological object and transmit a time sequence of images of the biological object; a camera located to receive the transmitted time sequence of images of the biological object from the microscope and transmit data representing the time sequence of images; a processor, coupled to receive the data representing the time sequence of images from the camera; the processor including an image processing algorithm that quantitatively determines mechanical deformation of the biological object from the time sequence of images, wherein the mechanical deformation is due to an interaction of the biological object with a plurality of molecules, wherein the image processing algorithm determines mechanical deformation by using an edge detector to determine an edge of the biological object and a time sequence tracker for tracking a plurality of movements of the edge by measuring positions of the edge at different times, and wherein the image processing algorithm also determines molecular binding kinetic values of the biological object as proportional to the mechanical deformation of the biological object.
 2. The system of claim 1 wherein the time sequence tracker for tracking the plurality of movements of the edge of the biological object comprises a differential image intensity analysis algorithm.
 3. A method for quantitative detection and analysis of the binding of molecules with molecular receptors on the surfaces of a biological object, the method comprising: operating a microscope and a sensor to capture a time sequence of images of a biological object; operating a processor to quantitatively determine the mechanical deformation of the biological object from the time sequence of images due to interaction of the biological object with a plurality of molecules, wherein determining the mechanical deformation comprises operating an edge detector for determining an edge of the biological object and tracking a plurality of movements of the edge by measuring positions of the edge at different times; and operating the processor to determine molecular binding kinetic values of the biological object as proportional to the mechanical deformation of the biological object and to determine binding kinetic constants, including ka, kd and KD from the time sequence of images.
 4. The system of claim 3 wherein tracking the plurality of movements of the edge of the biological object further comprises running a differential image intensity analysis algorithm on the processor.
 5. The system of claim 3 wherein the biological object is a cell and the image processor operates to measure the plurality of movements of the edge of the biological object by: detecting a cell edge from an optical image of the cell; selecting a region of interest (ROI) at a location of the cell edge; dividing the ROI into two equal halves, wherein one half is inside of the cell, and wherein the second half falls outside of the cell; and determining the edge movement by computing (A−B)/(A+B), where A and B are the intensities of the two halves of the image.
 6. A system for quantitative detection and analysis of binding of neuron stimulator molecules with molecular receptors on the surfaces of a neuron, the system comprising: a microscope positioned to image the neuron and transmit a time sequence of images of the biological object; an image sensor located to receive transmitted images of the neuron from the microscope and transmit data representing the time sequence of images; a processor coupled to receive the data representing the time sequence of images from the image sensor, the processor including an image processing algorithm that quantitatively determines mechanical deformation of the neuron from the time sequence of images, wherein the mechanical deformation is due to interaction of the neuron with a plurality of neuron stimulator molecules, wherein the image processing algorithm determines mechanical deformation by using an edge detector to determine an edge of the biological object and a time sequence tracker for tracking a plurality of movements of the edge by measuring positions of the edge at different times, and wherein the image processing algorithm also determines molecular binding kinetic values of the neuron as proportional to the mechanical deformation of the neuron.
 7. A method for quantitative detection and analysis of the binding of neuron stimulator molecules with molecular receptors on the surfaces of a neuron, the method comprising: operating a microscope and a sensor to capture a time sequence of images of a biological object; operating a processor to quantitatively determine mechanical deformation of the neuron from the time sequence of images due to interaction of the neuron with a plurality of neuron stimulator molecules, wherein determining the mechanical deformation comprises operating an edge detector for determining an edge of the biological object and tracking a plurality of movements of the edge by measuring positions of the edge at different times; and operating the processor to determine molecular binding kinetic values of the neuron as proportional to the mechanical deformation of the neuron. 